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Abstract 

We present a number models describing the sequential deposition of a mixture of particles 
whose size distribution is determined by the power-law p(x) ~ ai a_1 , x < I . We explicitly 
obtain the scaling function in the case of random sequential adsorption (RSA) and show 
that the pattern created in the long time limit becomes scale invariant. This pattern can 
be described by an unique exponent, the fractal dimension. In addition, we introduce an 
external tuning parameter (3 to describe the correlated sequential deposition of a mixture 
of particles where the degree of correlation is determined by (3, while (3 = corresponds 
to random sequential deposition of mixture. We show that the fractal dimension of the 
resulting pattern increases as (3 increases and reaches a constant non-zero value in the limit 
(3 — > oo when the pattern becomes perfectly ordered or non-random fractals. 
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1 Introduction 



The formation of stochastic fractals is an active field of research both theoretically and 
experimentally. Yet the mechanism by which nature creates fractals and the relationship 
between the degree of order and the fractal dimension is poorly understood. The history of 
describing natural objects by geometry is as old as the science itself. However, traditionally 
Euclidean lines, squares, rectangles, circles, spheres etc. have been the basis of our intuitive 
understanding of the geometry of almost all natural objects. But, nature is not restricted 
to Euclidean space. Instead, most of the natural objects we see around us are so complex 
in shape that conventional Euclidean space is not sufficient to describe them. It appears to 
be essential to invoke the concept of fractal geometry to characterize such complex objects 
quantitatively. It further enables us to search for symmetry and order even in disordered, 
complex systems |2|]. The importance of the discovery of fractals can hardly be exagger- 
ated. Yet, there is no neat and complete definition of a fractal. Instead one associates a 
fractal with a shape made of parts similar to the whole in some way. It is typically quantified 
by a non-integer exponent called the fractal dimension that can uniquely characterize the 
structure. This definition immediately confirms the existence of scale invariance, that is, ob- 
jects look the same on different scales of observation. To understand fractals, their physical 
origin and how they appear in nature we need to be able to model them theoretically. This 
forms part of our motivation of the present work. 

The simplest way to construct a fractal is to deterministically repeat a given operation. 
The construction of a classical Cantor set is a simple text book example of such a fractal. 
It is created by fragmenting a line into n equal pieces and removing n — m of the parts 
created and repeating the process with the m remaining pieces |]]. This process is repeated 
ad infinitum. However, this construction is rather too artificial as it differs in two ways from 
the fractals that occur in nature. It does not have any kinetics and it does not have any 
randomness. 

In this work, we introduce a stochastic process allowing a number of intrinsic tuning 
parameters which may be considered as a natural kinetic counterpart of the classical Cantor 
construction and which is a potential candidate in order to understand the essential governing 
rule of creating complex objects. These intrinsic tuning parameters are used to determine 
the degree of randomness and the rate at which a given operation is repeated to create a 
fractal. The interval are chosen for breaking stochastically and once an interval is chosen 
the cuts are placed randomly on the interval while the degree of randomness is determined 
by the precise choice of deposition kernels. Thus starting with an infinitely long interval, 
in the long time limit what remains in the long time limit are an infinite number of points 
scattered over the intervals. The properties of these points create a set that appears to be 
statistically self-similar and is characterized by a fractal dimension. 

The construction of stochastic fractals we consider is not at all pedagogical. One im- 
mediate and potential application of the stochastic fractal is the sequential deposition of a 
mixture of particles with continuous distribution of sizes. However, the model we consider in 
this work mimics the configuration when objects once inserted are clamped in their spatial 
positions for which non-equilibrium configurations are generated. The RSA processes have 
been found to describe many experimental systems, namely, the adhesion of proteins and 
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colloidal particles to uniform surfaces [[3], ffl, the reaction of various polymer chain systems 
such as methyl vinyl ketone and many fields in chemistry and physics. Although the 
process is conceptually simple, understanding its kinetics analytically is a challenging prob- 
lem (see an excellent review article ||). The deposition of particles of definite sizes in one 
dimension has been solved exactly and analytically in both continuous and discrete cases. 
The continuous version of this model is known as the random car parking problem. Recently, 
the deposition of a mixture of particles, such as deposition of two particles of different sizes 
has been considered H, H, |9| ■ From the experimental point of view, the deposition kinetics of 
well-defined mixtures of different types of particles is of great importance. It is well known 
that if objects are of definite size and once deposited are clamped in their positions, non- 
equilibrium configurations are created with a strong non-ergodic and non-Markovian nature 
||10|| . In this case, the resulting system does not produce a scale invariant pattern. Instead, it 
reaches a jamming limit when it is impossible to place further objects without overlapping. 
Hence a unique number (the jamming limit or the coverage) is sufficient to characterize the 
resulting pattern created in the long time. If a mixture contains a continuous distributions 
of sizes and is deposited sequentially, such that once deposited it is clamped to its position, 
then it is equivalent to our stochastic kinetic Cantor construction ||12|| . Hence the deposition 
of a mixture of particles will create a scale invariant pattern that does not reach the jamming 
limit since there always exists a particle if there is an uncovered space during the process. 

Interesting questions arises from the present work: (i) what is the role of the fractal 
dimension during pattern formation? (ii) Is there any relation between fractal dimension 
and degree of order? (iii) what are the relevant parameters to tune the degree of order 
and what are their physical meaning. The present work is an attempt to answer these 
questions. In section one, we present the general equation to describe sequential deposition 
of an arbitrary number of particles. In section two, we present model I describing random 
sequential deposition of a mixture of particles at a rate determine by the size of the empty 
space and exponent 7. We obtained explicit scaling solution and shown the resulting pattern 
in the long time limit is scale invariant both in time and space. In section three, we consider 
correlated sequential deposition of a mixture of particles (model II). We further obtain 
the fractal dimension of the resulting pattern for model II and compared with the fractal 
dimension for model I. Finally, we present a generalised version of both model to describe 
correlated sequential deposition of an arbitrary number of particles where, the degree of 
correlation is controlled by a parameter (3 and obtained fractal dimension as a function of 
(3. In the final section, we summarise the results discuss various point in order to reach a 
conclusion. 



2 Sequential deposition of arbitrary number of parti- 
cles 

The connection between the 1-d model of car-parking and the fragmentation processes was 



first emphasized by Ziff [ITTJ . The fragmentation process can be thought of as the deposition 



of points on a line whose position depends on the kinetic rule defined by the choice of kernels. 
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Let ip(x, t) be the concentration of empty or uncovered intervals of length x at time t. Then 
the rate equation for this concentration obeys the integro-differential equation 



-ip{x,t) I J[ p(x j )F(x n ^ m+1 ,..,x n \x 1 ,..,x n - m )S(x-J2x i )J[dx i + Tn 

j=l i=l i=l 

/n— 1 n— 1 n—m 
ip(y,t)F(x,x n - m+h ..,x n -i\x h ..,x n - m )5(y - x - J2 x i) II dxidy II P( x i) 

i=l i=l j=l 

(1) 

Here, n = 3,4,..., and m = 2,3,...,etc. so that (n — m) = 1, 2, 3, ...,etc is the number of 
particles which are deposited at each time step. Notice that the m value in particular put 
a strong constraint on the depositing particles. That is, if (n — m) = i and m < (i + 1) 
then the precise m will determine how many of the depositing particles must deposit next 
each other so that they produce only m new empty spaces. The term p(xi) determine the 
size of the depositing particles at each time step and -F(x n _ m+1 , .., x n \xi, .., x„_ m ) determine 
the rate and the rules with which Xi, .., x n - m particles are to be deposited to create m new 
empty spaces at each time step. The first term on the right represents the destruction of 
spaces of size x and the second term represents their creation from bigger spaces. Equation 
(1) can also describe the process of breaking an interval into n pieces and throwing away 
(n — m) pieces to create a stochastic fractal. 



dtp(x, t) 
dt 



3 Model I 

We first choose to consider the deposition of one particle at each time step (m = 2), when 
the size of the depositing particles follows the power law form 

p(x) = ax ' 1 for x < I. (2) 

This form of the parking distribution implies the deposition of a mixture of particles with 
a continuous distribution of sizes. We further choose the rate with which particles are 
deposited to be 

F(X2, X — Xi — X2\X\) = X 1 . (3) 

This choice of kernel implies that any point in the empty spaces are equally likely to be 
chosen for deposition by a particle whose size is determined by the choice of p(x)i However, 
the empty space where the particle is deposited is determined by the exponent 7. The rate 
equation then becomes 

-i>{x,t) + 2 / ^{y,t)y~<{y-x) a dy for x < I. (4) 



dt a + 1 

Notice that setting a = describes the random deposition of zero sized particles on a line 



at a rate x 7+1 i.e. the fragmentation process [13]. It is well known that in one dimension the 



random parking of cars of length '1' is highly non-ergodic in the sense that the whole space 
is not visited by the depositing particles. However, in the case of the random deposition 
of a mixture of particles with a power law form of distribution sizes, the system retains an 
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ergodic nature. In fact, the size of the particles to be deposited depends intrinsically on the 
size of the available empty spaces. Consequently, the size of the deposited particles becomes 
point like in the long time limit and there is always an available particle to deposit whatever 
the size of the empty space in the system. The second term on the right hand side of Eq. 
(4) reveals that as a increases the non-local character becomes more and more prominent, 
so that it becomes increasingly difficult to solve the equation. Notice that the exponent 7 
does not play any role in determining the degree of non-locality. 
We define the moment of the empty size distribution function as 

M(q,t) = t ip(x,t)x q dx (5) 
Jo 

We now multiply on both sides of the rate equation by x q and integrate over x thus obtaining 
a rate equation for the moments 

dM^t) = .!■(, + l)T(g + 1) 1 , 

This equation can be solved to find the solution for the qth moment, the solution taking 
the form of a generalised hypergeometric function with (a + 1) numerator and denominator 
parameters. In order to understand the various physical aspects of the problem within the 
simplest possible way, we first consider the case when a = 1. In this case the rate equation 
becomes 

8ib(r t) r a+y+l A 

-VM) + 2 / iP(y,tW(y-xrdy, (7) 



dt a + 1 

This class of models includes the one that is considered in |12|, when 7 = 0. Iterating 
the above equation to get all the higher derivatives and substituting in the Taylor series 
expansion of M(q, t) about t — yield 

M(q,t) = l q 2 F 2 (i^, i±^±l; i±I i±i; -tr»% (8) 
vy ' ; V7 + 2' 7 + 2 '7 + 2'7 + 2' /' v; 



where, 2-^2(0, b; c, d; x) is the generalised hypergeometric function |15], J6| and a = (—3 + 
\/T7)/2 = 0.5615288. The asymptotic expansion of the generalised hypergeometriv function 
immediately reveals that the moments show the following power-law behaviour 



q~~a 



M(q,t)~tT+*. (9) 

Notice that the exponent of the asymptotic expression for the moment is linear in q which 
reveals the existence of simple scaling. 



3.1 Explicit scaling solution 

In this section we attempt to show that in the long time limit the empty size distribution 
shows power-law behaviour. Linear power-law behaviour of the moments reveals that the 
system reaches to scaling or self-similar behaviour does not depend on the initial conditions 
for which one can invoke universality. However, a closer look at the rate equation further 
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reveals that only one of the two parameters has independent dimension i.e. x and t z have 
the same dimension where the exponent z is yet to be determined. That is, the dimension of 
ip must be expressible in terms of the independent parameter t or x alone. We can therefore 
define the dimensionless quantities as 

f = xT* or C = tx z (10) 

and 

«rfi) = *^>, *( te *) = ^. (id 

If scaling theory is obeyed a plot of against £ or a plot of against £ should fall on 
a single curve for any initial distribution. This reveal also that a self-similar solution in time 
and space exists. However, we find it convenient to consider the following scaling ansatz 

i/>(x,t) =x- e <5>{tx z ) (12) 

for which the pattern that develops in the long time limit is self-similar in space. The 
exponent 9 can be obtained from the rate equation for the moments using the condition for 
moment is time independent to give 9 = 1 + a. We now substitute this ansatz into the rate 
equation to obtain 

t * = — — — . (13) 

Demanding the scaling solution to exist we find that z = 7 + 2. The equation that we need 
to solve to find the scaling solution is 

7+2-9 1 7+2-9 2 f°° e_ 2 1 f°° 9+1 

C $'(0 + -C $(C) = — / v ~<+ 2 ®(v)dv r^C^ 2 / ^ -> +2 ®(v)dv (14) 

In order to eliminate the integral we differentiate this equation twice with respect to ( to 
reduce it to the third-order differential equation 

c v»(f) + c((3-^) + £K(C) + ((i-^)(i-^) + 5l±^cHc) 

We can rescale the equation to obtain 

<V (0 + C ((3 _ g+i, c w (0 + ((1 _ _g_ )(1 _ £+1 } + 37 + s-»Cw 

4 vs; 2V V 7 + 2 ; 2^ 7 + 2 A 7 + 2 ; 2(7 + 2) 2^ vs; 

+ f(l — )(1-^±1)- - 4 - Wfl=0. (16) 

V v 7 + 2 A 7 + 2 ; (7 + 2) 2 ^ vs; v ; 

The solution of equation (15) is given by the generalised hypergeometric function 

*(C) = 2 F 2 (l + ^" 2g - 1 i- ^ + 2g + 1 .i--g- l-?±i.-^ (17) 
lU 2^1 + 2(7 + 2) ' 2(7 + 2) ' 7 + 2' 7 + 2' ir { } 
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This is the explicit and exact scaling solution from which one can find the large-£ behaviour, 

$(C)~e-^ 2 . (18) 
The two scaling functions are related through 

(j)(xt^) = (xt^2)~ 9 $((xt^2) 7+2 )). (19) 
Hence we obtain the scaling function for large £ = £ 7+2 as 

<K0 ~ r (1+a) e-^ (20) 



One can recover the solution obtained in using different and indirect method from 
this general solution by setting 7 = 0. It is possible to obtain scaling solution for higher 
value of a but as a increases the numerator and denominator parameter of the generalised 
hypergeometric function becomes increasingly complicated. However, the knowledge of the 
z and 9 and details survey reveals that it is possible to write the scaling solution in the large 
£ limit for general a as 

<K0~r (1+a(Q)) e-^^, (21) 
where, a(a) is the solution of the following equation for q 

T(q + l)T(a + l) 1 , s 

2— — — - ' - = (22) 

T(q + a + 2) a + 1 y J 

Therefore, we can write the empty size distribution for the long time limit as 

i){x,t) ~x- (1+a(a)) $(C). (23) 

That is, we can choose scales tpo(x) = x~( 1+a( - a ^ depending on spatial variable for the empty 
size distribution function and to(x) = x~( a+7+1 ^ for the temporal variable. Therefore, in the 
new scale the properties of the empty size distribution function can be expressed in terms 
of one variable i.e., 

ip{x) ~ x- {1+a{a)) . (24) 
It implies that ip/ipQ and t/to are self-similar coordinates. 

3.2 Statistically self-similar pattern formation 

The existence of scaling shows that the pattern created in the long time limit becomes scale 
free i.e. the whole can be obtained from the part by a suitable change in scale. Essentially, 
this implies that we can invoke the idea of a fractal dimension: a dimension that uniquely 
determines the geometry of the object. We now use the usual box counting method to 
determine the fractal dimension. We define size of the segments to be 

i = S4~r*, (25) 

M(0,t) ' v ; 
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and we count the number of such segments needed to cover the whole set of points to 
determine the fractal dimension. In the limit 5 — > 0, we find that the number of segments 
< N(5) > requires to cover the set created by the Eq. (4) scales as 

< N(5) >~<T^ (Q) . (26) 

Where, Df(a) is the real and positive root of the polynomial equation in q obtained from 
the Eq. (6). Consequently, Df is found to obey 

r(,D / (a) + a + 2) a + l' 

Thus, a single scaling exponent Df(a) completely characterizes the structures of the objects 
which is reminiscent of the jamming limit. Note that fractal dimension Df(a) does not 
depend on the exponent 7 so it is independent of the rate at which particles are deposited, 
provided 7 > — (a + 1). In Table 2 we give a spectrum of fractal dimensions for different 
values of a with (3 = 0. This shows that as a increases the fractal dimension decreases. 
Later, we also attempt to give a physical interpretation of the exponent a. 



4 Model II 

We shall now consider another model that follow the same parking distribution (i.e. p(z) ~ 
but different deposition rate. The deposition rate of this model is 

F(x, y\z) = xyz (28) 

This particular choice of the deposition rate implies that all the points along the chosen 
empty space are not equally likely to be deposited. Although all the empty spaces compete 
on equal footing to be chosen where particles can be deposited. That is, the rate depends 
on the size of the deposited particles as well as on the size of the two smaller empty spaces 
created due to deposition. Substituting this into the equation Eq. (1) with n = 3 and m = 2 
we obtain the following rate equation 

0r{ - r - t] ■ v '" > -ip( x ,t) + 2 [ l x(y - x) a+2 i}(y,t)dy for x < I. (29) 



dt (a + 3)(a + 4) 

Notice that a = — 1 describe the fragmentation process with the fragmentation kernel 



F(x,y) = xy [ 19| . Also notice that the intensity of non-local character is higher than 
the previous model for which we find it is increasingly difficult to find scaling solution. Nev- 
ertheless, for our purpose it is enough to know that scaling exist. Substituting the definition 
of the moment into the above equation yields, 

dM(n,t) r r(n + 2)r(a + 5) 



dt 



T(n + a + 5) 



- 1 



M(n + a + A,t). (30) 



We find that the asymptotic behaviour of the moment can provide some of the interesting 
features of the system. Hence, from now on we are only interested in finding the fractal 
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dimension of the system that can uniquely characterize the structure. Following a similar 
procedure as we have done before we find that the number of segments needed to cover the 
whole set created, scales as 

< N(S) >~ 5- D ^ a \ (31) 
where, Df(a) as before can be obtained to satisfy 



ri 


[D f (a) + 2] 


T(a + 5) 


ri 


( K D f (a) + a + 5) 


l 



Here, the Df(a) is the fractal dimension of the pattern formed by this model. It is clear 
that the order of the polynomial equation is determined by the a value. Therefore the 
fractal dimension Df(a) is different for different a values. Comparing this model with the 
previous model we find that fractal dimension for this model is always higher than that for 
the previous model for each corresponding a value. Also, in both cases fractal dimensions 
appear to decrease monotonically as a increases. This is a feature that we shall discuss 
further later. We intend to determine if there exists any relation between the degree of 
order in the pattern and the fractal dimension. In order to do this we consider a further 
generalization of the two models we discussed. 



5 Model III 



We now turn to the more general model in which, at each time step, more than one particle 
will attempt to be deposited. We choose the parking distribution and deposition rates which 
are of the same functional form i.e., we choose the particle size distribution for deposition 
to take the form 

p(xi) = g(xi) for x < I, (33) 

and the deposition rate 

n 

F {Xn—m+1) • n •En\X\ ) ■ ■ i X n — m ) W dip^i)' (^4) 

i=n— m+1 

Substituting (29) and (30) into the equation (1) we obtain the following rate equation 
3ip(x, t) 

— — = -F n (x)if)(x,t) +mg(x) J ip{y 1 t)F Tl ^i{y-x)dy for x < I, (35) 
where, the functions {F n (x),n = 2,3...} are defined by 

/n 
5(x - ^2xi)Y[g(xi)dxi. (36) 
i=l 

Eq. (35) is equivalent to the dynamic system of breaking an interval into n pieces and 
throwing away (n — m) of them at each time step j^TJ . It is straightforward to show that 



F n (x) = [ g(y)F n -i(x -y)dy, (37) 
Jo 
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for n > 3 and 

F 2 {x) = [ X g(y)(x-y)dy. (38) 
Jo 

We further specify our model by choosing g(y) = y 13 , with j3 treated as an external tunable 
parameter. When g(y) = y 13 we can obtain the function {F n (x),n = 2,3, ..} from equations 
(32) and (33) as 

F ( X ) = TO + ^l" n(/?+l)-l /ggN 

It is interesting to note that if we set (5 = 0, m = 2 and n = a + 2 in equation (35) we get the 
same rate equation for the empty size distribution as described by the model I (equation (4)) 
with 7 = 0. It reveals that we can give an alternative interpretation of the model I. That 
is, at each time step an interval is broken into a + 2 random pieces and a of them removed 
from the system. Alternatively, we can say that the exponent a determines the number of 
particles to be deposited at each time step on an interval. However, for a > 1 the to value 
put a strong constraints on the depositing particles. That is, the process describe sequential 
deposition of a particles consecutively as if they were a single particle. Similarly, if we set 
f3 — 1, to — 2 and n = we get the same rate equation as described by the model II. This 
reveals that at each time step 9 ^ are deposited consecutively. These two feature helps us 
to understand the physical role played by the parking distribution exponent a. That is, as a 
increases the length of the depositing particles on the average increases with respect to that 
of corresponding lower a value. In order for further support we concentrate on the fractal 
dimension of the resulting pattern. 

Since the moments of the empty space distribution function can cheracterise the fragment- 
ing systems more easily than the empty spaces distribution function itself, we now consider 
the behaviour of the moment only. For g(y) = y 13 , the time evolution of the moments can 
be obtained using (5) 

This equation can be solved for the moments M(q,t), the solution again taking the form 
of a generalised hypergeometric function with (n — l)(/3 + 1) numerator and denominator 
parameters. We now consider the scaling behaviour of these models that define the scaling 
exponent 9 and z, and give the long time dependence of the moment as M(q,t) ~ t z ^ e ^ q ^ 1 \ 
We can immediately find z for all to, n and a, because, in the long time limit, the moments 
behave as 

M(q,t) ~ A(q)r b{q) (41) 
Substituting this into the rate equation for the moment yields a difference equation 

b(q + n((3+l) -1) =b(q) + 1. (42) 

Iterating this and using the q value for which the moment becomes time independent, we 
find 

%) = i^Tiyrr («) 
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This gives 

z = nfcg+Vi - (44) 

and q(P, m, n) related to 9 by 9 = Df + 1 and can be obtained from (36) and satisfies 

mT(q + (3 + l) T(P + 1) , . 

T(q + n((3 + l)) r(n(/3 + l))" 1 ' 

In particular, this model can be described as a correlated sequential deposition of particles 
on a substrate, where, the degree of correlation is determined by the exponent (3. 

For power-law form of parking distribution function this model describes the deposition 
of (n — m) = 1,2,3, ... etc. particles and creating m — 2,3, ...etc. number of new empty 
spaces respectively. If (n — m) = p and m < (p + 1) then the m value determines how many 
of them deposit consecutively. This deposition phenomenon can equivalently be interpreted 
as cutting an interval into n pieces and removing [n — m) of the parts created and repeating 
the process with the remaining m pieces thus resembling the concept of classic Cantor set. 
The fact that the size and position of the particles to be removed are chosen stochastically, 
unlike in the Cantor set where there are no kinetics, is irrelevant because we are considering 
the scaling regime, t — > oo. Since this process is repeated ad infinitum, it forms stochastic 
fractals with dimension between < Df < 1. As before, we use the box counting method 
by defining a characteristic length 5 so that we can count the number of segments required 
to cover the set when 5 — > which determines the properties of the resulting set and scales 
as, 

< N(5) >~ 8~ D f. (46) 

Thus we find that the Hausdorff-Besicovitch dimension in this case is equal to q*. We still 
have z given by equation (41), but now the value of q* (and hence 9 and Df) are non-trivial. 
This class of fractals includes that considered in |12|], with n = 3, m = 2 and (3 = 0. In the 
limit (3 — > oo we can use Stirling's formula and the following asymptotic formula 

T(az + b) ~ V2^e- az (az) az+b ~^, (47) 

in equation (42) to find that q*{= Df) — > ln(m)/ln(n). This coincides with the fractal 
dimension of the classic Cantor set. That is in the limit (3 — > oo the standard deviation of 
the size of the pieces created tends to zero. This is easily verified by calculating the mean 
and standard deviation from F{xi}. This particular finding implies that in this limit F n (x) 
behaves approximately as 

5(x\ — xs)dxi / 5(x\ — X2)dx2--- \ b{x\ — [x — ^ Xi))dx n -i. (48) 

i=l 

Substituting this into equation (32) we obtain 

d^fo*) = -—ip( x , t) + m(nx) x ip(nx, t). (49) 
at n 

This describes a model that splits an interval into n equal pieces and keeps only m of them. 
Substituting the definition of the moment Eq. (5) into Eq. (46) we obtain 

dM(q,t) _ r m li 
dt L n <?+i n . 



M(q + X,t). (50) 
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The solution of the Eq. (4.44) is 



M(q,t)~t- q -^, (51) 

and immediately reveals that Df — as in the classic Cantor set with kinetic exponent 
■r. In the limit (3 + 1 — > we can analyse Eq. (42) to show that q* (and Df) tends to zero 
like 7(771, n)(/3 + 1) where 7(m, n) = n(m — l)/(n — m). Consequently, we see that for all m 
and n there is a spectrum of fractal dimensions between (3 — > — 1 when Df — > and /5 — > 00 
when — > ln(m)/ ln(n). This is a very striking result. It implies that in the limit (3 — > 00 
particles are only deposited in the center of the empty space and produce strictly self-similar 
patterns. 

6 Summary and discussion 

In this work we presented a number interesting results relating to random and correlated 
sequential deposition of a mixture of particles of finite sizes. We find that the pattern 
created is statistically scale invariant. We also attempt to show the relationship between 
the fractal dimension and the degree of order in the resulting pattern created in the long 
time limit. In table 2 we present some values of the fractal dimension as a function of m 
and n for (3 = using model III. The corresponding fractal dimension for the Cantor set 
are given in parentheses. In table 1 we give the fractal dimension for m = 2 and n = 3 for 
some different values of (3. Table 1, and a more numerical survey confirm that the fractal 
dimension increases monotonically as (3 increases. Moreover, fractal dimension appears to 
decrease as n increases for a given m and vice versa, provided (n — m) > 0. In table 3, 
we compare the fractal dimensions for different values of a for model I and II. This table 
and further detail numerical survey confirms that for the same a value, the model II creates 
pattern that has higher fractal dimension than that for model I. We further notice that the 
first row of table 2 for fractal dimension with m = 2 are exactly the same as in table 3 
with (3 = and a = 1,2,3. Similarly, the fractal dimension obtained from model III for 
(3 = 1, m = 2 and n = would be the same as obtained from model II for corresponding 
= 1, 2, etc.. These, further shows that as a increases the length of the deposited 
particles on the average increases with respect to that of a corresponding average length 
for a lower a value. Which appear to be consistent with our details survey that reveals 
that fractal dimension monotonically decreases as a increases and as a — > 00 the fractal 
dimension Df —>■ 0. These observations immediately confirms that the exponent a does not 
play any role in creating ordered pattern since fractal dimension does not reach to constant 
non-zero value as a — > 00. 

In a recent letter |T7|] , Brilliantov et. al studied the random sequential adsorption of 
a mixture of particles with a continuous distribution of sizes determined by the power-law 
form (equation (2)). They reported that the pattern created in the long time becomes more 
and more ordered as a increases and in two dimensions it reaches the Apollonian packing in 
the limit a — ■> 00 when the depositing particles are finite mixture of disk. It is important to 
notice that the one dimensional analogue of this model is deposition of a mixture of rods. 
Evidently, one expect the fractal dimension to coincide with the classic Cantor set when 
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the pattern becomes perfectly ordered pattern. In Ref. [I7| an exact expression for fractal 
dimension also given for general d dimension. It is also clear that the Apollonian packing 
is a non-random or strictly self-similar fractal like Sierpinsky gasket difference lies only in 
the geometry of the depositing particles. The observation that the pattern becomes more 
and more ordered should be true for any dimension for any geometry including d — 1. That 
is, in order to support the result that in the limit a — > oo the pattern becomes more and 
more ordered the fractal dimension must reach to constant non-zero value in that limit. In 
particular, as a — > oo the fractal dimension must reach to the value In 2/ In 3 which is one 
dimensional analogue of Apollonian packing and Sierpinsky gasket both. In one dimension 
we can solve the model exactly which corresponds to our model I. The exact enumeration 
of fractal dimension and details numerical survey reveals that Df—>0asa^oo instead 
of reaching a constant non-zero value. Therefore, it contradicts the analysis we give in this 
work with result reported in JT7j. More recently, deposition of a finite mixture of rectangles 
in two dimensional substrate is studied in Ref. ||20|| . Although, in this work particles are 



only allowed to deposit one of the four corner, yet, it retains the generic feature of deposition 
phenomena of a definite mixture and in particular it is very close to the model we consider 



in this work. The work in ref. [2D| is stochastic counterpart of the deterministic Sierpinsky 
carpet or Cantor gasket which is strictly self-similar. In this deterministic case, the 
initiator is a square and the generator subdivides at each step into b 2 equal square of which 
p of them are removed according to a fixed rule. After an infinite number of iterations 
the resulting set can be seen as a generalization of the Cantor set to two dimensions that 
has the fractal dimension Df — ln(b 2 —p)/hxb. However, in the case of its stochastic 
counterpart we have shown that the system does not reach to simple scaling instead the 
system shows multiscaling. A result which we believe to be true for deposition of a finite 
mixture particles of any geometry in more than one dimension. That is, the pattern created 
in the long time limit has global scaling exponent Df and local scaling exponent known as 
f(q), where q is the Holder exponent. That is, pattern can be divided into a subset that 
scales with different fractal dimension, a phenomena called multifractality. However, in the 
case of strictly self-similar pattern the system reaches to simple scaling behaviour. The study 
further revealed that the global exponent or the fractal dimension of random fractal is always 
lower than its corresponding strictly self-similar counter part. Therefore, we conclude that 
fractal dimension must increases with increasing order and reaches to maximum value when 
the pattern is in perfect order. In this work, we show that the exponent a does not play 
any role in creating an ordered pattern. Instead, it implies that the length of the deposited 
particle at each time step increases on average as a value increases. 

In order to create an ordered pattern we reveal that one needs to choose F(x, y\z) ~ (xy) 13 
and p(z) ~ z@ . This model for (3 ^ describes that at each event an empty space is chosen 
randomly. However, once this decision has been made, particles are more likely to deposit 
at the center of the empty space than on either side of it. Of course the tendency to deposit 
at the center increases as /3 increases. In fact, our analysis further supports that the fractal 
dimension increases with the degree of increasing order and reaches its maximum value (a 
non-zero constant) in the perfectly ordered pattern, as it does in the classic Cantor set or in 
the Seirpinsky gasket. Krapivsky et al. reported in [12] that the dimension of the random 
fractal is always smaller than its deterministic counterpart. If we look at the situation for 
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the stochastic Cantor set, we find that in order to get the deterministic classic set (j^f), 
a alone does not play any role in creating any ordered pattern. Instead, one has to choose 
F(x,y\z) to be the same power law form as for parking distribution p(z). In this case only 
the a value determines the degree of tendency to place the particles in the center of the 
empty space and in the limit a — > oo particles are always placed exactly at the center of the 
empty space. We can quantify the increasing regularity of the resulting pattern created in 
the long time limit by introducing the concept of entropy production that characterizes the 
degree of order as 

S = -Y,P(Ck)hgp(C k ). (52) 

Since in the f3 — > oo limit there is only one definite configuration we have p{Ck) = 1 that 
contributes to the entropy. In fact we can define 

p(C k ) = = D \ , (53) 
mm/ Inn 

where, each Df corresponds to one definite configuration. We find that S increases as 
(3 decreases towards zero. That is, as (5 decreases the number of possible configurations 
increases. In this work, we generalize the conventional RSA where the position of a particle 
to be deposited in the empty space is chosen randomly and that the degree of randomness 
by the position dependent deposition rate. That is, the position where the particle is to be 
placed is chosen by the size of the empty space being destroyed and by the size of the empty 
spaces created on either side. As a prospect of future work one can choose p(z) = 8(x — 1) and 
F(x,y\z) to be position dependent. Studying RSA with these choices the system obviously 
will reach a jamming limit, but how it varies with the degree of order can be of greater 
physical interest. 
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Table 1: The fractal dimension Df of stochastic fractals for m = 2 and n = 3 for increasing 
/3 values. 
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Table 2: The fractal dimension Df of the stochastic fractals with (3 — 0; the corresponding 
dimensionality for the Cantor set {(3 = oo) is given in the parentheses. 
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0.5616 
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0.5183 
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0.466542 
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0.33405 


0.429121 
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0.30784 


0.400614 
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0.288505 
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7 


0.27351 
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Table 3: The fractal Dimension Df for (3 — 0, 1 for different values of a 
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